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Abstract: In this talk I review some recent applications of NRQCD. I first discuss the unquenched 
NRQCD lattice extractions of the strong coupling constant, paying particular attention to the re- 
cent advances in reducing systematic errors. I then discuss the progress made in testing the 
NRQCD/factorization formalism for onium production. In particular, I gather the evidence, or lack 
thereof, for the universality of the production matrix elements. I also discuss the interesting question 
of polarized production, which is a crucial test of the formalism. I address the viability of this theory 
for the J/tp system as well as what needs to be done before we can reach any definitive conclusions. 



The study of heavy quarks is interesting from 
a phcnomcnological as well as formal standpoint. 
Indeed, the utility of the experiments designed 
to study the CKM sector of the standard model 
is bounded by our theoretical understanding of 
states containing heavy quarks. For theorists 
this presents a great challenge given that, in gen- 
eral, it is quite difficult to make predictions for 
strongly interacting particles from first princi- 
ples. However, we do have one powerful tool at 
our disposal which saves us from despondency, 
namely effective field theories. These theories 
represent limits of QCD which make approximate 
symmetries manifest, thus greatly enhancing our 
predictive power. Furthermore, since the effec- 
tive theory reproduces QCD in a well defined 
limit, it is possible to calculate corrections in a 
systematic fashion. 

The effective theory I have been charged to 
review is non-relativistic QCD. As the moniker 
implies, this is a theory which approximates full 
QCD when applied to a bound state containing 
more than one heavy quark. The quarks neces- 
sarily have small velocities in the limit jtlq 3> 
A-qcd, due to asymptotic freedom. The theory 
is written down as a simultaneous expansion in 
as{2rriQ) and the relative velocity in the cen- 
ter of mass frame of the quarks, v. Indeed, if 
we wish, we may calculate corrections to arbi- 



trary order in powers of v and a s , at the price of 
introducing unknown, yet universal parameters. 
Thus, the relevant question is not "is NRQCD 
correct?" , but rather, "how well do the v and a s 
expansions converge for the system of interest?" 

In this review I will discuss recent advances 
in the application of NRQCD. In particular, I will 
concentrate on applications for which we have 
new data. Thus, due to space limitations, I will 
not discuss the interesting topic of inclusive onia 
decays, nor will I discuss recent theoretical progress 
on threshold top quark production ffl]. Instead, I 
will focus on the utilization of NRQCD to mea- 
sure the strong coupling, and the predictions for 
onia production. 

1. The Effective Field Theory 

The effective Lagrangian is constructed by en- 
suring that at the matching scale, ttiq, the effec- 
tive field theory reproduces on shell Greens func- 
tions of QCD to some fixed order in a(2mg) and 
v. The operators in the effective field theory are 
then classified by their "effective v dimension" |s| . 
That is, to each operator we may assign a power 
of v which is determined by some velocity power 
counting rules. For our purposes, we will be 
interested in the fermion anti-fermion sector of 
the theory where the lowest order Lagrangian is 



Heavy Flavors 8, Southampton, UK, 1999 



Ira Rothstein 



given by 



f D 2 

L = / dtd 3 x^ f (iD + — )*(i,x) + 
7 2m 



x - y 



(1.1) 



where V' and x are the fermion and anti-fermion 
two spinors respectively, L g and L q are the un- 
scathed dimension four operators of the QCD La- 
grangian for the gluons and light quarks respec- 
tively. Cy is unity at tree level but is corrected at 
higher order. The coefficients of the other opera- 
tors in the heavy quark sector are fixed by repa- 
rameterization invariance [Q. Terms of higher 
order in v are easily calculated by expanding the 
full theory on-shell matrix elements in the v. In 
the one fermion sector this leads to the usual 
terms encountered in non-relativistic quantum 
mechanics, such as the magnetic moment and 
Darwin interactions. There is one technical de- 
tail which I should point out at this point. When 
matching beyond leading order in the strong cou- 
pling, the above Lagrangian must be modified if 
one wishes to preserve manifest power counting. 
In particular, if one wishes to regulates integrals 
using dimensional regularization, then the non- 
potential interactions must be dipole expanded^, 
rc| . That is to say, we must make the replacement 



A^t, x) A^t, 0) + x • dA^t, 0) + 



(1.2) 



This dipole expansion incorporates the fact that 
the typical gluon momentum is of order mv 2 whereas 
the size of the state is order mv. On a calcu- 
lational level this leads to gluonic interactions 
which do not transfer three momentum, which in 
turn ensures manifest power counting. The situ- 
ation is actually slightly more complicated that 
this, because there are also gluons which have 
typical momentum of order mv |?J . This leads to 
further modifications of the Lagrangian, which I 
will not have space to discuss in this review || |^| . 

2. Lattice Extractions of a s 

As far a lattice calculations are concerned, the 
motivation for effective held theories is different 
than the one discussed in the introduction. In 



particular, by removing the heavy quark mass 
from the theory it is possible to work with lat- 
tice spacings a, which are larger than the Comp- 
ton wavelength of the heavy quark (ma) > 1. 
This alleviates the problem of the quark liter- 
ally falling through the cracks. Indeed, the use 
of coarser lattices allows for precision measure- 
ments which would otherwise be presently un- 
manageable. 

Perhaps the most interesting use of lattice 
NRQCD is for the purpose of extracting the value 
of the strong coupling constant a s . This low en- 
ergy extraction is interesting, not only because 
it allows us to compare to results obtained from 
perturbative QCD, but also because any dispar- 
ity between the low energy and high energy ex- 
tractions, after renormalization group running, 
could imply the existence of new physics. Thus, 
a precise low energy extraction is of great scien- 
tific interest. Below we will discuss the results 
of the two collaborations, "NRQCD" fuj and 
"SESAM"||, which used NRQCD to treat the 
heavy quarks. These are the first lattice extrac- 
tions of a s performed in the unquenched approx- 
imation. While the NRQCD collaborations re- 
sults have been around for several years now, the 
important complimentary results of SESAM col- 
laboration are more recent. The two collabora- 
tions use different techniques so that their agree- 
ment would provide strong corroboration of the 
results for this important measurement. 

The strategy behind the extraction is to first 
define a short distance coupling constant a s (a) 
via some lattice observable. Both collaborations 
utilized the so-called plaquette coupling [jll) , which 
is related to small Wilson loops via 

-logW M = —a p (3Al/a) + 

(1 - a p (1.1879 + (0.0249, 0.070)71/)) .(2.1) 

The first number in parenthesis is for Wilson 
fermion and the second is for staggered fermions, 
and 



-(ReTr(pexp -ig s f A-alx 



.(2.2) 



W m , n is the Wilson around a rectangular path 
of size (ma) x (na). This observable has the fur- 
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ther advantage that its leading non-perturbative 
correction, 



SW„ 



-7ra 4 (mn) 2 
36 



(a s G 2 ). 



(2.3) 



is believed to be small due the anomalously small 
value of the gluon condensate (a s G 2 ) ~ 0.05 GeV 4 . 
Note there are no higher order perturbative cor- 
rection to the relation |2.1| , as a matter of defini- 
tion. This just shuffles higher order corrections 
into the relation between ajj-g and a p . The scale 



3.41/a is the BLM scale 13 determined by cal- 
culating the n f dependent piece of the two loop 
perturbative expansion of the Wilson loop. The 
lattice spacing is then determined by measuring 
the the spin independent quarkonia splittings. 
These splittings are chosen since they are much 
less sensitive to tuning of the bare heavy quark 
mass. Once the value of a p (3Al/a) has been ex- 
tracted, it can be related to otjjg via 



— ap + C 2 a P 

TT 



(2.4) 

The scale e 5 ^ 6 Q is chosen to absorb the BLM 
piece of the first order coefficient. The coefficient 
C2 is only known in the pure gauge theory E^], 
in which it takes the value Ci — 0.96. Using 
this relation its possible to determine a value of 
ajj S (Mz) by running with the three loop beta 
function and taking into account the relevant quark 
thresholds. 

There are several sources of errors, and the 
SESAM and NRQCD collaborations incur them 
in different fashions, as I shall now discuss. As 
usual with lattice calculations there is the issue 
of cut-off artifacts. The NRQCD collaboration 
removed all 0(a 2 ) error from the heavy quark 
action by using an improved action. They did 
not improve the gluonic action to this order, but 
made a perturbative estimate of these correc- 
tions. They then checked that their spin split- 
tings are independent of the lattice spacing (in 
the quenched approximation) by calculating with 
several different spacings in the range, 0.05 — 
0.15 fm, and found no variations. The use of 
staggered fermions also introduces errors on the 
order of 0(a 2 ). They estimated the net effect of 
lattice artifact errors to be at the level of 0.2%. 



The SESAM collaboration used dynamical Wil- 
son fermions which incur lattice spacing errors 
at linear order in a. They worked at fixed lat- 
tice spacing and were therefore unable to per- 
form a scaling analysis. They quote a larger er- 
ror, due to discretization, than did the NRQCD 
collaboration, namely %5. This larger error can 
be attributed to the difference between the light 
fermion actions. 

In addition, there are corrections due to trun- 
cating the Lagrangian at some fixed order in v. 
Both collaborations include 0(v 2 ) corrections to 
the effective action. 0(v 4 ) corrections should be 
negligibly small in the T system since the 0(v 2 ) 
corrections were found to shift the mass splitting 
by 10%. There is also an error incurred by trun- 
cation the perturbative expansion of the Wilson 
coefficients for these higher dimension operators. 
The coefficients of these operators arc tree level 
"tadpole improved". An estimate for this error 
was made by the SESAM collaboration by vary- 
ing the tadpole improvement prescription. They 
found that this error is of the same order as the 
relativistic errors, and estimated these combined 
errors to be at 1% level. Whereas the NRQCD 
collaboration estimated the relativistic errors to 
be at the 0.2% and the errors due to tadpole pre- 
scription to be 0.5%. 

The use of dynamical light quarks was a cru- 
cial step in making these lattice predictions trust- 
worthy. However, it is difficult, at this time, 
to perform calculations with physical masses for 
the light quarks. Therefore, the calculations are 
performed with unphysical light quark masses, 
and the results are extrapolated to the physical 
value. In it was shown that the level split- 
tings should grow linearly with the light quark 
mass. The SESAM collaboration calculated the 
splittings for several different light quark masses, 
between the radial as well as orbital excitations. 
The results for the IS — IP splitting are shown 
in figure 1 and seem to fit a linear relationship 
quite well. The result gives us confidence that 
the extrapolation to physical light quark masses 
is being performed correctly. Both collaborations 
estimate the error due to use of unphysical sea 
quark masses to be at the 1% level. 

Both collaborations performed their unquenched 
calculations with n/ = 2, and extrapolated to 
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Figure 1: SESAM result for the 1 P 1 - 3 Si splitting 
as a function of the light quark mass. 

n/ = 3. Indeed, the NRQCD collaboration found 
that fixing the lattice spacing using different T 
splittings led to results which varied by 3 a when 
rif is taken to be zero. This discrepancy essen- 
tially disappears when rif — 3. So it seems rea- 
sonable to expect that these new unquenched cal- 
culations have the light quark effects well under 
control. We should point out that neither collab- 
oration took into account SU (3) violation , but 
these effects should be extremely small, given the 
size of the bound state. 

A larger error arises due to the uncertainty 
in the, yet to be calculated, nj dependent piece 
in the relation [0| The SESAM group varied 
the coefficient C2 in [2.4| , between 1 and -I and 
found a 2%— 3% variation in their result while the 
NRQCD collaboration quote a 1.9% error due to 
the uncertainty in C2 . The SESAM collaboration 
quotes the results 

( 5 ) I \ 

f 0.1118 (10)(12)(5) x - T splitting 

\ 0.1124 (13)(12)(15) T' - T splitting 1 ^ j 

where \ is the spin average of the P wave states. 
The first error is statistical while the last two are 
from the NRQCD truncation and uncertainty in 
the sea quark mass dependence. The errors do 
not include any guess at the size of discretization 
error stemming from the light quarks, which the 
authors estimate to be at the 5% level. Nor does 
this error include the uncertainties due to our 
ignorance of the n / dependent piece of C2 , which 
the authors estimate to be at the 2% — 3% level. 



The NRQCD collaboration quotes 

( 5 ) I \ 
a MS^ = 

f 0.1174 (15)(19) x-T splitting 

\ 0.H73 (21) (18) T' — T splitting 1 ' ' 

The errors are due to lattice artifacts and 
perturbative truncation errors, respectively. No- 
tice that the two collaborations differ by 3a. The 
SESAM collaboration attributes this to light quark 
discretization errors. 

We see that these calculations seem to have 
the errors well under control. My only true con- 
cern, is the issue of the convergence of the pertur- 
bative expansion at these scales. Given the pre- 
cision of the measurement, a calculation of the 
complete C2 coefficient would certainly increase 
the confidence level. 

3. Onium Production 

In Q it was pointed out that by combining per- 
turbative factorization with NRQCD it is possi- 
ble to make "rigorous" 1 predictions for onium 
production. A general production process may 
be written as 

da = d<? l+ ^ QQ[n]+x (0 I Ol I 0). (3.1) 

a 

Here da i+ ^QQ^ +x is the short distance cross 
section for a reaction involving two partons i and 
j, in the initial state, and two heavy quarks, in 
a final state labeled by n, plus X. This part of 
the process is short distance dominated and com- 
pletely calculable in perturbation theory, mod- 
ulo the possible structure functions, in the initial 
state and may be considered as "matching co- 
efficient" . The long distance part of the process 
involved the hadronization of the heavy quarks in 
the state n into the hadron of choice H . Indeed, 
the matrix element which is written as 

(O|o^|o} = (o|^ r »'x|^if + x) 

x 

(# + x| x t r>|o>. (3.2) 

The tensor T n operates in color as well as spin 
space and also contains possible derivatives. This 

lr rhe calculation have a level of rigor equivalent to 
those in semi-inclusive hadro-production where one relies 
on perturbative factorization in the physical region. 
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tensor also determines the order, in v of the ma- 
trix element. The size of the matrix element is 
fixed by determining the perturbations necessary 
to give a non-vanishing result for the time or- 
dered product (selection rules). The matrix ele- 
ments are taken between states of the effective 
theory, which are given by the eigenstates of the 
dipolc expanded version of 1.1. Thus, at lead- 



ing order there is no overlap between operators 
with quantum numbers (usually labeled spectro- 
scopically 2S+1 ij(l,8)) which differ from those 
of the state under consideration. Which is to say 
if H has the quantum numbers 2S+1 Lj then at 
leading order in v color singlet operators give 

(0 | Of^'+'L'j,) | 0) oc Sss'Sl^Sjj,, (3.3) 

while all color octet operators give zero. Higher 
order contributions can be included by inserting 
higher multipole moment interactions into the 
time ordered product. For instance, the matrix 
element (0 | J J? ( 3 Si) | 0), would scale as v , 
since we need two El insertions at the the ampli- 
tude level, each costing a factor of v. The first in- 
sertion neutralizes the color, but it also changes 
L. So a second insertion is needed to bring us 
back to the S wave state. This is illustrated in 
figure 1. 




Figure 2: Gluon fragmentation into a Si state. 
The quark anti-quark pair are created at a short dis- 
tance scale in a relative octet state. They then prop- 
agate into the J/ip state by emitting two soft El 
gluons. 

For the color singlet matrix elements, we may 
use factorization, which is typically correct up to 



v 4 corrections, 2 to write 

(o | of | o) = (o | ^r n x | h)(h I xtr-v | o). 

(3.4) 

We may then interpret (H \ x^T n ij} | 0) as the 
"wave function" at the origin ip(0), or derivatives 
thereof, in a potential model (with some factor- 
ization scale dependence). This parameter is the 
same parameter which determines the annihila- 
tion rate and thus its value is easily extracted 
from the data. 

Note that if we were to drop all color octet 
contributions, then we would end up with the 
old fashioned, color singlet model. In this model 
the production rate is calculated by weighting 
the partonic rate to produce two heavy quarks 
with zero relative velocity with the a potential 
wave function at the origin | ip(0) | 2 - Here I 
wish to emphasize that this method of calcula- 
tion is indeed a model, in that there is no limit 
in which it reproduces the full QCD calculation. 
On the other hand, despite what might often be 
said, the NRQCD prediction is NOT a model. 
In the limit where the quark mass is much larger 
than the QCD scale the theory reproduces full 
QCD. Thus, those who call the NRQCD predic- 
tion, the "color octet" model, have been either 
gravely mislead or are using a definition of the 
term model with which I am not familiar. 

The power of the NRQCD/factorization for- 
malism for production lies in the universality of 
the matrix elements (0 | Of 0). So that to 
make a prediction for a given process we must 
first extract the appropriate matrix elements from 
another process. This universality is a conse- 
quence of the factorization of the short distance 
production process from the relatively long dis- 
tance hadronization process. At large pt pro- 
duction is dominated by fragmentation processes 
|l7|] . In this case, the cross section is written as 
the probability to create a nearly on-shell parton 

2 , times 



with momentum q such that qo ^> \ 
the probability of the parton to fragment into the 
hadron of interest. This latter probability may 
be written in terms of an unphysical fragmenta- 
tion function Df 1 (z) which gives the probability 



2 This is true unless the L of the state created is dif- 
ferent from the L of the operator under consideration, in 
which case the correction can be 0(v 2 ). 
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of forming a hadron H with momentum fraction 
z from an initial parton labeled by i. This fac- 
torization is analogous to case of factorization 
for semi- inclusive light hadron productionjl6j . In 
the case of onium production the fragmentation 
function itself factorizes into a product of the 
probability for an energetic parton with invariant 
mass on the order of 4mQ to produce two heavy 
quarks in a state 2S '+ 1 £^ 1 ' 8 ) times the probabil- 
ity of this heavy quark state to hadronize into 
the onium state of interest. Thus, as opposed to 
the case of light hadrons, the fragmentation func- 
tion can be written in terms of a perturbative ex- 
pansion in as(2mQ) and an unknown constant (s) 
(our matrix elements from above). Writing the 
cross section in this way makes it simple to re- 
sum large logs of the form \ogpr l m Q (i-e. per- 
form the Alterelli-Parisi running). In this case, 
the factorization holds up to corrections of order 
Am 2 Q /p T . 

For smaller values of px fragmentation no 
longer dominates and one must consider the fu- 
sion of partons into the final state. This cross 
section is written as the probability to create a 
heavy quark anti-quark pair, in a given configura- 
tion, times the hadronic matrix element. In this 
regime higher twist effects are expected to scale 
as powers of A/ \fp\ + m-QV 71 3 , though this has 
not been fully explored to date (n has not been 
fixed), since the quark pair can be strongly in- 
fluenced by the jet in the forward direction as pt 
gets smaller. Thus, our confidence in our calcu- 
lation dwindles as pr is reduced. We should thus 
keep in mind that when we test the formalism we 
are testing more than just the convergence of the 
a s and v expansions. We are also testing factor- 
ization. Thus, at least for production, and to a 
lesser degree for decays, we should really say that 
we are testing the NRQCD/factorization formal- 
ism. 

3.1 The Universality of the Matrix Ele- 
ments 

A first crucial test of the NRQCD/factorization 
formalism consists of checking the universality 
of the matrix elements. There are many places 



3 The first power correction would most likely scale as 



from which we may extract the relevant matrix 
elements. The first extractions were done us- 
ing the Tevatron data. Let us, for the moment 
concentrate on the case of the J/ip, where we 
presently have most of our data. As discussed 
above, at large pr gluon fragmentation domi- 
nates and gives a leading order contribution pro- 
portional to (0 | Og( 3 S{) | 0), as shown in 
figure 2, while at intermediate values of of pr, 
the rate becomes sensitive to (0 | Og^So) | 0) 
and (0 | Og /,/ '( 3 Po) I 0). Actually at lower values 
of pt the rate is sensitive to a linear combination 
of these last two matrix elements. This combi- 
nation is usually defined as 

M m = {0 >/*£so)) + A (0 ^ ( 3p o)) , (3. 5) 

m Q 

k actually varies with px- The extraction of 
is complicated by the fact that its value is very 
sensitive to the small x gluon distribution which 
makes it difficult to extract reliably. Varying 
parton distribution function can change the ex- 
tracted values by a factor of four. Therefore, I 
will concentrate on extractions of (0$^ ( 3 Si)} . 

Several groups have extracted this matrix el- 
ement Jy| |cj H |D for the J/r/> as well as 
the ip 1 . I have collected the various extracted 
values from these references in tables 1 and 2. 
The extractions differed in several respects. Ref- 
erences [^0], [21], 22 included a Monte Carlo es- 
timation for initial state radiation, which we ex- 
pect to become more important for smaller val- 
ues of pr- Of these three references though, only 

correctly accounted for the Alterelli-Parisi 
evolution, which is seen to be important for this 
particular matrix element. Noting the difference 
between J22j and [^0[ |2l|, we can see that the 
effect of resummation is only important, as we 
would expect, for small pt- Only |l9| made any 
estimates for the scale dependence of the results, 
which from the table can be seen to be substan- 
tial. The extractions also differ in the way in 
which they treat the interpolation between the 
fragmentation at high pt and the direct produc- 
tion at low pt- I think that it would reason- 
able to say that the dominant source of error 
will come from higher order corrections, as ev- 
idenced by the large scale dependence found by 
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J/1> 


CTEQ2L 


CTEQ4L 


CTEQ4M 


GRVLO 


GRVHO 


MRS(R2) 


MRSDO 


111 














0.66 ±0.21 


1- 




1.06 ±0.14™ 




1.12 ±0.141°;" 




1.40±0.22iJ;^ 






0.33 ±0.05 








0.34 ±0.04 




0.21 ±0.05 


y 






0.27 ±0.04 












0.96 ±0.15 








0.92 ±0.11 




0.68 ±0.16 



Table 1: Extractions of (Og /l/, ( 3 Si)> in units of 1CT 2 GeV 3 . The first error is statistical, while the second, 
when listed, is due scale dependence. 





CTEQ2L 


CTEQ4L 


CTEQ4M 


GRVLO 


GRVHO 


MRS(R2) 


MRSDO 


















0.46 ±0.10 




m 




0.44 ± 0.081^^4 




o.46 ± omtojl 




0.56±Q.ll±g;g 






20| 


0.14 ±0.03 








0.13 ±0.02 




0.11 ±0.03 



Table 2: Extractions of {Of ( 3 Si)) in units of 1(T 2 GeV 3 . The first error is statistical, while the second, 
when listed, is due scale dependence. 



|19j . It would be interesting to see the extent 
to which this scale dependence is reduced by in- 
cluding higher order effects. However, given the 
present state of affairs, I don't think I would be 
overly conservative if I were to say that there is 
a factor of two uncertainty in the octet matrix 
element. 

At LEP, we might hope that we can get a 
better handle on the octet matrix elements since 
the theoretical calculation is perhaps under bet- 
ter control. In this case we need not worry about 
factorization scale dependence in parton distri- 
bution functions or initial state gluon radiation. 
However, there is a complication do to final state 
soft gluon emission. Indeed, J/ip production at 
LEP is a beautiful example of color coherence. 
Moreover, as opposed to hadronic collision, soft 
gluon effects are completely calculable in closed 
form in lepton initiated processes, as I will know 
discuss. 

If one calculates the leading order differential 
cross section for octet production |2^, Q one 
finds schematically 

^-(Z - J/i> + X)^a 2 s log(Mf/M^)/z + , 

(3.6) 

where z — 2Ej i^/Mz and the terms left off are 
less singular in the z — > limit. We see that the 
differential rate is dominated by an infrared en- 
hancement coming from the small z region. This 
enhancement leads to large double logs in the 



total rate. Given that cvlog(M|/Mj/^,) ~ 1, a 
leading order prediction necessitates a resumma- 
tion of these logs. Such a resummation can not 
be accomplished by standard ladder resumma- 
tion techniques, due to the color coherence of 
the soft gluon radiation [^5|. However, by im- 
posing angular ordering it is possible to write 
down an integral equation for the resummed frag- 
mentation function p6[ . The resummed differ- 
ential rate calculated in |27j] is shown in figure 
3. Notice that the total rate is dominated by 
the octet at small z. The peak at small z is 
an example of what is called the " hump-backed 
distribution" that arises in calculations of jet- 
multiplicities |2l| . The prediction is not valid 
at smaller values of z due to the complete break- 
down of the perturbative expansion, as well the 
the breakdown of a saddle point approximation 
used in (27). Nonetheless, the region well below 
the peak contributes negligibly to the total rate. 
The data from LEP includes feed-down from ex- 
cited states, so that the rate is not proportional 
to (Og ), but instead is proportional to a combi- 
nation of matrix elements. This is the drawback 
of this study. The rate is in fact proportional to 
the combination 

m 

BR{ip(m) -> J/tjj + X), (3.7) 
where ip(m) are excited states of the J/ip. The 
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authors of IE 71 found 



<0 8 J/,A ( 3 Si)) = (0.019±0.005 stat ±0.010 tWa ) GeV' J 

(3.8) 

The theoretical error quoted is conservative, and 
includes the uncertainty due to relativistic cor- 
rections, subleading logs and factorization scale 
dependence. If one compares this with an anal- 
ogous extraction from the Tevatron data, which 
has errors on the order of 100%, one finds agree- 
ment. Before moving on I should point out that 
the peak is almost purely color octet. Thus the 
data seems to support the existence of this chan- 
nel with a matrix element of order that found at 
the Tevatron. 



T 

es i 




Figure 3: Results for dT/dz at LEP taken from 
pTj. The dashed line shows the result without the 
resummation. 

Before concluding our discussion of univer- 
sality, I should address the so-called "HERA- 
anomaly" . It has been pointed out that at HERA 
in the large z 4 region, NRQCD predicts a sharp 
rise in the inelastic J/ ip direct photo-production 
cross-section . This behavior is due to the 

fact that the 1 Sq and 3 Pj configurations can be 
produced via a t-channel gluon at lowest order in 
a s . A similar statement can be made of the 3 Si 
configuration for resolved photo-production. No 
such peak near z « 1 is seen in the data leading 
to the aforementioned "anomaly" . However, as 
is now well appreciated, as z approaches 1, the 
theory for the spectrum becomes intractable for 
several reasons. Firstly, the endpoint region is 



sensitive to initial state intrinsic transverse mo- 
mentum (equivalently initial state radiation [^3| ). 
Indeed, the authors of [M found that by mod- 
eling the intrinsic transverse momentum with a 
Gaussian distribution with (kx) ~ 0.7 GeV they 
were able to fit the data. Secondly, as z ap- 
proaches its endpoint (z > .75) the spectrum be- 
comes sensitive to the long distance hadroniza- 
tion process. This manifests itself as a break- 
down in the NRQCD expansion In j3§ 
it was shown that the rate may be written in the 
following convoluted form 



da 

dz 



p-> o 

T,min 



5 s(l — z)x — 



dP% J dx J dy+ S{x, z, P|) F H [n] (y+) 
Af 2 (l-z)+P 2 

(3.9) 



M 2 (l-zf+P 2 



z{l-z) 



-y+ 



Ej/^/E-y in the protons rest frame 



The lower cut on px is necessary to eliminate the 
diffractive contribution. The function Fn[n](y+) 
is a leading twist distribution function defined in 
which physically accounts for the momentum 
carried by the non-perturbative gluons. For any 
non-zero Pt we see that the expansion parameter 
close to z = 1 is — z) ~ w 2 /(l — z). Thus, 

the NRQCD expansion breaks down for 1 — z. ~ 
0(v 2 ), because higher-order terms in v 2 grow 
more and more rapidly as z, — > 1. Consequently, 
the NRQCD factorization approach makes no pre- 
diction in the endpoint region and the discrep- 
ancy between leading order predictions and the 
data in this region does not allow us to draw any 
conclusion on the relevance of color-octet con- 
tributions to photo-production. If one averages 
the z-distribution over a sufficiently large region 
containing the endpoint, the octet mechanisms 
contribute significantly to this average. However, 
the characteristic shape information is then lost 
and one has to deal with the more difficult and 
uncertain question of whether the absolute mag- 
nitude of the cross section requires the presence 
of octet contributions, and whether their mag- 
nitude is consistent with other production pro- 
cesses. One might hope to learn more about the 
octet matrix element from studying the smaller 
z. However, there are large uncertainties in the 
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color singlet contribution 35 , 36) which make the 
prospect of learning much from this region rather 
bleak. 

3.2 The question of polarization 

NRQCD predicts that at large transverse mo- 
mentum, when gluon fragmentation dominates, 
the onium should be dominantly transversely po- 
larized [[38]. This polarization arises as a con- 
sequence of the fact that the gluon is nearly on 
shell and thus is transversely polarized up to cor- 
rections of order Am 2 c /p T . This transverse po- 
larization is inherited by the onium, since the 
electric transitions preserve the spin of the heavy 
quarks. There are corrections coming from hard 
gluons which can flip the spin but these correc- 
tions turn out to be surprisingly small |59). There 
are also spin symmetry breaking magnetic tran- 
sitions which are suppressed by v 4 compared to 
the leading order transition. Thus, we have the 
rather robust prediction that ip production at 
large pr should be dominantly transversely po- 
larized and this polarization should increase with 
Pt- 

Presently there is data from CDF for both 
polarized J ftp and ip' production [[IT]. The data 
on the J/ip shows no signs of polarization. How- 
ever, this data includes feed-down from higher 
excited states, which we would expect would di- 
lute the polarization 5 . The ip' sample is perhaps 
more surprising. In this case we don't expect to 
have the dilution problem, and thus if our ap- 
proximations are good, we should see a rise in 
the fraction of transversely polarized ip n s as pr 
increases. The preliminary data does not seem to 
show this trend. Indeed, it seems to indicate a 
trend towards longitudinal polarization with in- 
creasing pt- Of course, the data still has large 
errors and thus it is perhaps as little too early to 
jump to any conclusions. However, the trend is 
rather disturbing. 

4. Conclusions and Outlook 

In this review, I have touched upon only two ap- 
plications of NRQCD. There are several other ap- 
plications which were discussed at this meeting 



f|. As far as using NRQCD for extracting a s , I 
think that the next important step will be to try 
to get a better handle on the perturbative rela- 
tion between cijjg and ap, which will be a rather 
Herculean task. However, given that we believe 
that the series is asymptotic, there will always 
be doubts as to how well the series is behaving, 
and given the accuracy of these calculations, one 
always worries. 

As far as the NRQCD/factorization formal- 
ism, as applied to the charmed system, is con- 
cerned we are presently in a state of uncertainty 6 . 
Let us review the successes and failures of NRQCD 
in this context. To begin with let's consider pre- 
dictions for cross-sections. The data at the Teva- 
tron for ip' seems to necessitate a color octet con- 
tribution [Q. The spectrum is well fit, once the 
octet contributions are included. This is not a 
terrible surprise since there are two octet con- 
tributions one which behaves as l/p T the other 
as l/pfr, nonetheless this is encouraging. The 
overall normalization of the prediction is uncer- 
tain at about a factor of two level. At LEP the 
theory is under better control, but the data is 
much more sparse. The theory predicts a large 
hump at small energies due to the octet contri- 
bution which seems to be supported by the data. 
A combined LEP analysis of the spectrum would 
certainly be welcomed. As far as the HERA data 
is concerned the data can be well fit without 
the contribution from the octet channel at next 
to leading order [j36[ The expected rise at 

large z from the octet, which is not seen, should 
not trouble us at this time, for reasons discussed 
above. Octet contributions to production in fixed 
target experiments 43, [iJL f45fl , as well as in B 
decays UJ] have been analyzed. These stud- 
ies also seem to yield matrix elements of order 
10 -3 GeV 3 . However, we would expect higher 
twist contributions in these reactions to be espe- 
cially important. 

As far as the qualitative size of the extracted 
matrix element is concerned, NRQCD seems to 
have it right. If we use the relation B.4 to ex- 
tract the singlet production matrix element from 
a decay process, we find that the ratio of octet to 



3 A quantitative study of this dilution would be helpful. 



°I have not discussed the T system due to space limita- 
tions and the fact that until recently the data was rather 
scant. 
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singlet matrix elements is of order 10 -3 . Given 
that we expect v 2 ~ 0.3 in the J ftp system and 
that according to NRQCD scaling rules the ratio 
of these matrix elements should scale as u 4 , we 
seem to have qualitative agreement. However, I 



would like to point out that the relation 3.4 itself 



has never been tested. Which is to say, we have 
yet to extract the size of the singlet production 
matrix element. Indeed, all extractions of octet 



matrix elements assume the relation 3.4 and then 



use the decay matrix element, either from a po- 
tential model, or from an NRQCD extraction. 
One could worry that the production singlet ma- 
trix element is actually smaller than expected. 
As far as I can tell, this can not be ruled out by 
the data, at least for the J /if) 7 . One easy way to 
check, would be to extract the singlet production 
matrix element at CLEO, where the singlet dom- 
inates, away from the endpoint . While I don't 
expect the singlet to be anomalously small, it is 
still important to extract the singlet production 
matrix element, since the perturbative series for 
decays seems to be very poorly behaved [ blj . 

For more quantitative results, we should stick 
to the large pt processes where we feel more con- 
fident that higher twist effects are small. This es- 
sentially leaves us with the Tevatron data and the 
LEP data. The Tevatron extraction is plagued 
by large factorization scale dependence, which 
could indicate a poorly behaved perturbative se- 
ries. While the next to leading order correction 
to the total cross-section has been calculated , 
a next to leading order calculation of the differen- 
tial cross-section da/dpx would allow us to study 
higher order corrections without worrying about 
higher twist contamination. The LEP extrac- 
tion, on the other hand, is limited by statistics. 
Thus, at least for now, testing the universality of 
the matrix elements, at an accuracy level beyond 
100% is difficult. 

Finally, there is the issue of polarization. In 
addition to the studies at the Tevatron, polar- 
ization has also been investigated in fixed tar- 
get experiments, where the theory is compati- 
ble with the experiment within errors, given the 



7 Here I am not including the data from fixed targets, 
since the higher twist effects could be large. Also, at 
HERA for z 0.5, it might be difficult, but given the 
uncertainties, not impossible. 



large uncertainty in the octet matrix elements. 
There have also been interesting theoretical stud- 
ies of J / ip polarization in B decays |32| , LEP J5(J 
and at HERA for photo-production |35| as well 
a lepto- production [ p3| , which could in the fu- 
ture provide useful information. Presently, we 
should be concerned with the trend in the Teva- 
tron data. If, once the statistics improves, the 
data continues to show no trend towards trans- 
verse polarization at larger px then we have a 
very intriguing puzzle on our hands. Where could 
the theory be going wrong? A naive guess would 
be the spin symmetry is badly violated. Which is 
to say that perhaps the magnetic transition op- 
erators are not as suppressed at we think. How- 
ever, while possible, this seems unlikely given the 
fact that it seems to work so well in the D meson 
system. If spin symmetry were badly violated, 
then we would expect the splitting between the 
D and the D* to be much larger than it is. Of 
course, the matrix element which determines this 
splitting is different than those which induce the 
spin flipping transitions in ip' production. So we 
can't draw any hard conclusions from smallness 
of the splitting. The next naive guess would be 
that there are large perturbative corrections to 
the fragmentation which lead flip the spin. But, 
as discussed above, these corrections have been 
calculated |43) , and are indeed small. Another 
possibility is that gluon fragmentation does not 
dominate at large pr- For instance, it could 
be that the production singlet matrix elements 
(Of ( 3 Si)) is 30 times larger than expected from 
relation ft.4[ The problem with this proposal is 
that it would mean that the LEP data is too 
small by an order of magnitude at large energies. 
Thus, there is no obvious way to depolarize the ip' 
at large pt, no less get them to be longitudinally 
polarized. I should stress, that the prediction 
for polarization is based on very basic assump- 
tions, which have been tested in other context 
many times. Those assumptions are: Factoriza- 
tion at large px, spin symmetry, and the stan- 
dard parton model assumptions which are rou- 
tinely tested in Drell-Yan processes. These as- 
sumptions are "derivable" from QCD, and have 
been tested repeatedly. So if the data persists 
we have a real puzzle on our hands. Much more 
light will be shed on this problem as more data 
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streams in. With more upsilons we may hope 



[13 



to gain information as well, since in that sys- 



tem we believe all of our expansions will be bet- 
ter behaved. However, we will have to go to 



[14 



larger values of pt to see the expected transverse 
polarization [59 1 . Nonetheless we still should still 



[15 



see the trend. 
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